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ABSTRACT 



We estimate the stochastic gravitational wave (GW) background signal from 
the field population of coalescing binary stellar mass black holes (BHs) through- 
out the Universe. This study is motivated by recent observations of BH-Wolf- 
Rayet star systems and by new estimates in the metallicity abundances of star 
forming galaxies that imply BH-BH systems are more common than previously 
assumed. Using recent analytical results of the inspiral-merger-ringdown wave- 
forms for coalescing binary BH systems, we estimate the resulting stochastic 
GW background signal. Assuming average quantities for the single source en- 
ergy emissions, we explore the parameter space of chirp mass and local rate 
density required for detection by advanced and third generation interferomet- 
ric GW detectors. For an average chirp mass of 8.7M Q , we find that detection 
through 3 years of cross-correlation by two advanced detectors will require a rate 
density, r > 0.5 Mpc" 3 Myr _1 . Combining data from multiple pairs of detectors 
can reduce this limit by up to 40%. Investigating the full parameter space we find 
that detection could be achieved at rates ro ~ 0.1 Mpc~ 3 Myr _1 for populations 
of coalescing binary BH systems with average chirp masses of ~ 15M Q which are 
predicted by recent studies of BH-Wolf-Rayet star systems. While this scenario is 
at the high end of theoretical estimates, cross-correlation of data by two Einstein 
Telescopes could detect this signal under the condition r > 10~ 3 Mpc _3 Myr _1 . 
Such a signal could potentially mask a primordial GW background signal of 
dimensionless energy density, fi GW ~ 10~ 10 , around the (1-500) Hz frequency 
range. 

Subject headings: gravitational waves - binaries: close - cosmology: miscella- 
neous 
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Introduction 



Coalescing systems of stellar mass binary black holes (BBHs) are among the most likely 



candidates for the fir st detection of gravitational waves (GWs) (IFlanagan fc Hughes! Il998 
Buonanno et al.ll2003l ). Their enormous predicted luminosities ~ 10 23 L^, would allow fu ture 



ground-based int erferometric GW dete ctors, such as Advanced LIGO (IHarry et al.ll20f0l ) and 
Advanced Virg o (lAcernese et al.l 120 09) or third generation instruments such as the Einstein 
Telescope (ET; iPunturo et al.l 120101 ). to probe these sources o ut to Gpc distances. In this 
paper, we are motivated by recent increased rate estimates (IBelczynski et al.l l2010bl ) to 
explore the possibility that a population of BBHs could form a detectable stochastic GW 
background (SGWB) signal for these instruments. 

SGWBs can result from the superpositio n of populations of unresolved primordial 
( iGrishchuk! 11974 ) or astrophysical sources (see iRegimbaul l201ll for a recent review). As- 
trophysical SGWB signals are important for at least two reasons. Firstly, they contain rich 
information on the global properties of source populations, such as their source rate evolu- 
tion, their mass ranges and their average energy emissions. Secondly, a dominant continuous 
astr ophysical background could mas k the relic SGWB signal from the very early Universe 



sec 



Maggiordl2000t iBuonannol 120031 . for reviews). 



Mer gers of binary neutron stars have been suggested as sources of potential ly detectable 
SGWBs ( IRegimbau fc de Freitas Pachecoll2006l ; IRegimbau fc Chauvineaull2007l) . Recent ob- 



servations of BH-Wolf-Rayet (WR) star systems (ICrowther et al 



2010) and new estimates 



in the metallicity abundances of star forming galaxies ( IPanter et al.l 120081 ) imply that the 
Galactic merger rate o f BBHs may be of a similar order to that of binary neutron stars 
( IBelczynski et al.ll2010bf ). Therefore, a population of more luminous BBHs could produce a 
dominant background signal. Our aim is to explore upper limits for a SGWB from coalescing 
BBHs over a range of rates and system masses. We investigate the constraints future ground 
based interferometric GW detectors will be able to place on the average properties of the 
BBH population. 

The paper is organized as follows. In Section 2 we discuss rate estimates of coalescing 
BBHs and then derive cosmic source rate evolution models for different star formation his- 
tories and minimal delay times. In Section 3, source ene rgy spectra for coa lescin g BBHs are 



obtain ed using the template gravitational waveforms of lAjith et al.l (120081 ) and lAjith et al. 
( 120091 ). We then calculate the BBH background in Section 4 and discuss the detection 



regimes, detectability and constraints on the parameter space of the predicted background 
in Section 5. Finally, in Section 6 we present our conclusions. 
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Rates of BBH coalescences 



The discovery of BH-WR star systems within close proximity has increased rate esti- 
mates of coalescing BBHs. There are presently two known systems: NGC300 X-l, which 
l ies at a distance of 1 .8 Mpc and is composed of a ~ 2OM BH and a WR star of ~ 26M© 
dCrowther et al-lboioh : IC10 X -l contains a BH of a m ass at least 23M and a ~ 35M Q WR 
star, and lies within 700 kpc (IPrestwich et al.l 120071 ). As WR stars are the progenitors of 
Type Ib/c supernovae, if such systems survive the supe rnova explosion, B BH systems will 
form and eventually coalesce within a timescale of Gyrs (IBulik et al.l 1201 ll ). 



Recent results from Sloan Digital Sky Surve y, have indicated t hat half of recent star 
formation involved galaxies with low metallicity (IPanter et al.ll2008l ). This has a profound 
effect on the coalescence rates of compact binaries containing BHs when one considers that 
NGC300 X-l and IC10 X-l were both formed in low metallicity environments. 

Survival of BBH systems is highly dependent on whether they can overcome two key 
obstacles in their evolution. Firstly, post natal supernova kicks can disrupt a significant pro- 
portion of systems. Secondly, orbital shrinkage during the common envelope phase when the 
larger star transfers mass to its smaller companion , can cause the stars to merge before they 
become compact objects. iBelczynski et al.l (l2010fi ) have shown through population-synthesis 
modeling that a lower metallicity environment can suppress these two effects. Firstly, obser- 
vational evidence suggests that l arger BH masses, which are produced at low metallic ity, are 
born with lower kick velocities (IMirabel fc Rodriguesl l2003t IBelczynski et al.l l2010al ). Sec- 
ondly, in lower metallicity environments, slower radial expansion occurs during the common 
envelope phase, thus increasing binary retention. The greater fraction of systems that can 
survive, in combination with a greater detection range from more massive and hence lumi- 
nous systems, h as increased the detectio n prospects of BBHs for ground-based interferometric 
GW detectors faelczvnski et alJboiObh . 



Previous estimates of the coalescence rate of BBHs formed through isolated binary evo- 
lution in the field have ranged over orders of magnitude, from 10~ 4 to 0.3 Mpc~ 3 M yr -1 with 
a realistic value of 5 x 10~ 3 Mpc~ 3 Myr _1 jKaloeera et al.ll2007t kbadie et allhoioh . The ef- 
fect of metallicity discussed above increases the realistic estimate to 3.1 x 10 -2 Mpc _3 Myr _1 
assuming a 50-50 mixture of solar and 10% solar metallicity and the mos t stringent evolu- 
tionary scenario with respect to system survival ( iBelczynski et al.ll2010bl ). Unless we state 
otherwise, this rate (denoted as r±) is adopted in our calculations. 

A higher estimate of 0.43 Mpc~ 3 Myr -1 (denoted as r 2 ) was obtained by assuming that 
all systems survive early merger during the common envelope stage. We note that r 2 leads 
to a detection rate for initial LIGO/Virgo of around 5 events per year. A recent population 
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synthesis study by iBulik et al.l (1201 if ), however, has demonstrated the formation of BBHs 
with high chirp masses (~ 15M Q ) from the two BH-WR systems, and has yielded a similarly 
high rate of 0.36Mpc~ 3 Myr _1 corresponding to 3.6 detections a year. They suggest that 
either currently employed searches are insensitive to higher mass BBH inspirals or that there 
is an additional aspect to the evolution of such systems that has not so far been considered. 
To take account of this uncertainty, we take r2 as a higher rate. 



2.1. The cosmic rate evolution model 



The cosmic coalescence rate can be extrapolated from the local rate density r by as- 
suming the rate tracks the star formation rate (SFR). Explicitly, the differential GW event 
rate in the redshift shell z to z + dz can be written as 



dV 

dR = roe(z)——dz, 
dz 



(1) 



with dV the cosmology dependent co-moving volume element given by 



dV 

dz 



Attc 



r(z) 2 
Wz) 



(2) 



where the Hubble parameter H(z) = Ho[Qa + Cl m (l + z) 3 ] 1 ^ 2 and r(z) is the comoving 
distance related to the luminosity distance by di = r(l + z). We use the param eters H = 
lOO/i ■ kms^Mpc" 1 with h = 0.7, fi m = 0.3 and fi A = 0.7 fcomatsu et al.ll2009h . 



Source rate density evolution is accounted for by the dimensionless evolution factor e(z), 
norma lized to unity in our local intergalactic neighbourhood. Following iRegimbau fc Hughes 
(120091 ). we define e(z) = p* jC (2:)//?* iC (0) normalized to unity at z = where 



P*A Z ) 



(1 + z f ) 



P{t d )dt d , 



(3) 



relates the SFR to the BBH coale scence rate. Here, p* is the SFR density in M yr _1 Mpc -3 , 
based on the parametric form of iHopkins fc Beacoml (120061 ). derived from recent measure- 
ments of the galaxy luminosity fun ction. To allow for u ncertainties in the SFR, we also 
consider a recent model described in I Wilkins et al.l ( 120081 ). obtained through measurements 
of the stellar mass density - this model gives a much lower rate for z > 1 . The factors z and 
zj represent the redshift values of BBHs merger and BBH system formation respectively. 
The delay time for BBHs, td, is given by the difference in lookback times between Zf and z 
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dzf 



(1 + z')H(z') 



(4) 



Here P(td) is the probability distribution of delay times. For this we take the form, 1/td, 
with a lower cutoff to at 100 Myr as suggested by current popu l ation synthesis studies 
(IBelczynski et al.l 120021 ; iPostnov Yungelsonl 120061 ; iDominik et al.l l201ll ) . To account for 
uncertainty in the lower cutoff and possible correlation between P(td) an d BBH chirp masses, 
we additionally consider a longer minimal delay time of t = 500 Myr (T. Bulik, private 
communication) . 

Figure 1. shows the cosmic rate evolution factor e(z) of BBH coalescences for different 
SFRs and minimal delay times. Although the models show some variation in peak z, we 
show in section 4 that the magnitude of the background is largely dependent on r rather 
than the form of e(z). 



3. The GW energy spectrum of coalescing BBHs 



The evolution of a coalescing BBH is traditionally divided into three phases: inspiral, 
merger and ringdown. While the early inspiral and ringdown phases can be approximated 
analytically by post-Newtonian expansion and perturbation theory, to model the late inspiral 
and merger requires a numerical solution of the Einstein equations. In the last few years 
breakthroughs in numerical relativity have enabled the inspiral-merger-ringdown evolution 
of B BH coalescences to be model ed with high accuracy for a broad space of parameters (see, 



e.g. 



Hannam 2009; Hinder 2010, for reviews). 



For the complete evolution history of coalescing BBHs, phenomenological waveforms 
can be constructed by frequency domain matching of p ost-Newtonian inspiral waveforms 



with coalescence wayeforms from numerical s imulations (lAjith et al.l 120071 ; iBuonanno et al. 



20071 ; iPan et al.ll2008l ; ISantamaria et al.ll2010l ). Such waveforms share a common feature, in 
that the Fourier amplitude is approximated to a leading order as a power-law function of 
frequency z/~ 7 / 6 for the inspiral phase, followed by v^ 2 ^ for the merger stage and a Lorentzian 
function around the quasi-normal mode ringdown frequency for the ringdown stage. 

For this study, we convert to a energy spectru m to account for th e individual source 
emissions. We choose two templat e mod els: one from lAjith et al.l (120081 ). for a n on-spinning 



case, a nd another from lAjith et al.l (120091 ) for BBHs with non-precessing spins. In lAjith et al. 



(120081 ) the Fourier amplitude is given by their equation (4.13), which we convert to an energy 
spectrum dE jdv. As we expect the inspiral spectrum to equal the post-Newtonian approx- 
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imation (see, e.g.. ICutler et al.lll993t iFinn fc Chernoflall993f ). the inspiral-merger-ringdown 
spectrum for a BBH with component masses m,\ and m<i is given by 



dE 
dv 



(Gtt) 2 / 3 M ( 



5/3 



-1/3 if 



V < V\ 



bJ\V 



UJ 2 



2 / 3 if^< 



V < v 2 



(5) 



l + (^) 2 . 



if V2 < V < V3 



Here M c is the chirp mass, M, 



5/3 



mim 2 mi + 7712] 



-1/3 



^ 1 and 0J2 



-4/3 



are constants chosen to make dEjdv continuous across V\ and The set of parameters 
{i>x, z/ 2 , cr, z/ 3 ) can be determined by the two physical parameters (the total mass M and the 
symmetric mass ratio rj) in terms of {arf + br] + c) / 7rM, with coefficients a, 6, c given in 
Table 1 of kiith et all j2008h . producing (404, 807, 237, 1153) Hz for a 1OM -1OM BBH 
(M c = 8.7M ). 



The waveform presented in lAjith et al.l (120091 ). includes spin effects through a single 
spin parameter x = (1 + 5)xi/2 + (1 — 5)x 2 /2, with 5 = (m.j — m 2 )/M and Xi = Si/mf. The 
parameter represents the spin angular momentum of the ith black hole. The corresponding 
Fourier amplitude includes a minor correction (related to x an d v) f° r non-spinning BBHs. 
We construct energy spectra for BBHs with non-precessing spins based on their equation (1). 

Figure 2 shows the GW energy spectra for a 1OM — 1OM BBH assuming: the non- 
spinning case; x = 0.85; x = and y = —0.8 5. The two extreme values for % are set by 
the numerical simulations of lAjith et al.l (120090 . corresponding to both binary components 
having maximal spins aligned or anti-aligned with the orbital angular momentum. The 
radiation efficiencies for these energy spectra are 6.7%, 9.74%, 5.15% and 4.28% respectively. 
We note that the radiated GW energy mainly depends on M c and Xi an d that the energy 
spectra for v < 100 Hz show little variation. We note that effect of orbital eccentricity is not 
considered in our derivation of energy spectrum. This has little e ffect on our results as the 
orbits of coalescing compact objects are expected to circularise (jPeterd Il964j ) before their 
GW signals reach the sensit ive frequency band of ground-based interferometric detectors 
( IBrown fc Zimmerman! |2010| ). 
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4. The stochastic GW background 

In this section we evaluate the spectral properties of the BBH background. It is cus- 
tomary to characterize the SGWB by the energy density parameter 

^Gw(^obs) = — ff GW , (6) 
p c d In z/ obs 

where pew is the GW energy density, z/ bs is the observed GW frequency and p c = 3Hq/8ttG 
is the present value of critical energy density required to close the Universe. For a SGWB of 
astrophysical origin, fi GW is related to the spectral energy density F v (in ergcm _2 Hz _1 s _1 ) 
by 

^Gw(^obs) = -r^F u (u ohs ), (7) 
where F v at the observed frequency z/ Q b s can be written as 

f u (is ohs ,z) — (z)dz, (8) 

Here, z sup = min(z max , f m ax/^obs — 1), with z max the maximum redshift of SFR model and z/ max 
the maximal emitting GW frequency. The differential GW event rate, dR/dz, is given by 
equation (1) and f v is the energy flux per unit frequency emitted by a source at a luminosity 
distance d^z) 

where dE/dv is the gravitational energy spectrum given by equation (5) and v is the fre- 
quency in the source frame which is related to the observed frequency by v = f bs(i + z). 

Figure 3 shows the energy density parameter Qqw of the BBH background corresponding 
to five models outlined in Table 1. The final model (e) approximates the common power- 
law behavior of models (a-d) before reaching a peak of ~ at 400-600 Hz, and uses a 
conservative estimate to account for effects of spin, SFR and delay time. This model is given 
by ^GW^obs) = 1-95 x 10~ n z/^s f° r 1 Hz< z/ obs < 400 Hz. We note that using the two 
extreme spectra (x = 0.85 and x — —0.85) in Figure 2 leads to a variation in the results of 
less than 40%. 
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We see in Figure 3 that effects of spin, SFR and delay time are insignificant for estimating 
the BBH background considering all other uncertainties in sources rates and systems masses. 
Combining the equations used to calculate fi GW , we find the BBH background is sensitive 

5 /3 

to r and average value of M c through the relation fi GW oc r M c ■ It is useful to express 
model (e) in the form: 



^Gw(^obs) ^ 4.2 x 10 



-10 



r 



3.1 x 10- 2 Mpc" 3 Myr~V V8-7M 



Mr 



5/3 



Vobs \ 

100 Hz 7 



2/3 



(10) 



allowing us investigate the parameter space (r , M c ) for a SGWB in the next section. We 
note that equation (10) is consistent with equation (10) of lPhinneyl (120011 ). Considering that 
the most sensitive frequency regime of the second and third generation detectors is within 
100-200 Hz, we further adopt a cutoff frequency of 400 Hz when applying this function. 



5. Issues of detectability 

A number of laser interferometric GW detectors have reached their design sensitivities 
and have been coordinating as a global array. These include the LIGcfl detectors based at 
Hanford (H) and Livingston (L) in USA, the Virgc§ (V) detector in Italy and the GEO60C§ 
detector in Germany. The LIGO and Virgo detectors are undergoing upgrades that will 
produce a order of magnitude improvement in sensitivity. Advanced LIGO is expected to 



1 http : //www . ligo . caltech . edu/ 
' http : //www . virgo . inf n . it/ 
http : / / www . geo600 . uni-hannover . de 



Model source spectrum SFR model Delay time, to (Myr) 



a 


non-spinning 


HB06 


100 


b 


x = o 


HB06 


100 


c 


non-spinning 


HB06 


500 


d 


non-spinning 


W08 


100 


c 


power law model 







Table 1: The 5 models and their parameters used to determine the energy density parameter 
foGW °f a BBH background. Model e, a power law approximation, is described in the text. 
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be operational by 2015 and an Advanced Virgo facility is expected to begin commissioning 
in 2011. Other plans for advanced detectors include LCGT0 (the Large-scale Cryogenic 
Gravitational- wave Telescope; denoted as C), an underground detector to be built by the 
TAMA group in Japan and LIGO-AnstraliJ (A), a 4km interferometric detector that would 
be based in Western Australia. Third generation detectors have also been proposed, such as 
ET, with a target sensitivity 100 times greater than current instruments. 

To evaluate the detectability of the BBH background, we consider the five advanced de- 
tectors introduced above (H,L,V,C,A). We take LIGQ- Austr alia to have the sa me sensitivity 



as Advanced LIGO and for ET adopt ET-B s ensitivity from iHild et al.l ([20081 ) as well as the 
latest ET-D sensitivity from lHild et al.l ( 201ll ). Design sensitivity curveqjfor these detectors 
are shown in Figure 4. 



5.1. Duty Cycle 

For a SGWB of astrophysical origin, in addition to the energy density parameter and 
characteristic frequency, another useful quantity is the duty cycle (DC). This is defined 
as the ratio of the typical duration of a single signal to the average time interval between 
successive events, e.g., 

DC = T(z)dR(z), (11) 
Jo 

where dR(z) is the differential event rate given by equation ([!]) and z max corresponds to the 
redshift limit of considered SFR models, f (z) is the average observed duration of GW signals 
generated by individual sources at redshift z, given at leading order by 



5c 5 



2567r 8 / 3 G 5 / 3 



[(l + z)M c }- 5/3 ^. 



12) 



Here, z/ m i n is the lower frequency bound of the detector, which is set by the low-frequency 
seismic "wall" for ground-based interferometric GW detectors. The present LIGO detector 
has v mm — 40 Hz, and this can be reduced to 10 Hz for advanced detectors and 1 Hz for ET 



sec 



Regimbau fc Hughes! 12009k for details). 



http : //gw. icrr .u-tokyo . ac . jp/lcgt/ 



5 http : //www . aigo . org . au/ 



3 Data are taken from |http: / / wwwcascina.virgo.infn.it / advirgo / 



bin/DocDB/ShowDocument?docid=2974, http://www.et-gw.eu/etsensitivities 



https://dcc.ligo.org;/ cgi- 
and iKuroda et alT|2010l ) . 
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In general, a GW background with a DC of unity or above is denned as continuous^. And 
non-continuous signals can be further categorized into popcorn noise (0.1 < DC < 1) and 
shot noise {DC < 0.1) type. As source rate evolution will increa se out to large cosmological 



volumes, some stu dies investigate how DC too increases with z fjCoward fc Regimbaul 12006 



Howell et al.ll201ll ). In this study we are concerned with the value of DC as seen at the 
detector and thus equation (jTTj) is interpreted as a total value. 

For a background signal produced by a source population with an average chirp mass 
(M c ) = 8.7 M Q , for rates r x and r 2 we find DC values of 0.01 and 0.2 for advanced detectors. 
Therefore the signal will most likely be of the shot noise category (or at most popcorn noise 
for the higher rate r 2 ). For ET type detectors the signal will be continuous with DC values 
of 5.8 and 80.6 for rates ri and r 2 respectively. 

Although for advanced detectors, the SGWB ca lculated in this pape r is no t continuous 
(Gaussian) even at the higher rate r 2 , we note that iDrasco fc Flanagan! (120031 ) have found 
the cross correlation method nearly optimal for a DC > 10 -3 over 1-year integration^!. 
In the following sections we therefore consider the cross correlation statistic to assess the 
detectability of the estimated BBH background. 



5.2. Cross correlation of detectors 



The optimum detection strategy for contin uous GW background s i gnals is cross-co rrelating 
the output of two neighbouring detectors (see. lAllen fc Romanolll999l ; lMaggiordl2000l ). This 
requires that the detectors are separated by less than one reduced wavelength, which is about 
100 km for frequencies around 500 Hz where ^gw(/) might peak. The detectors also need 
to be sufficiently well separated that their noise sources are largely uncorrelated. We note 
that although this may not be possible for ET, te chniques are in d evelopment to remove 
environmental noise and instrumental correlations (iFotopoulosl 120081 ) . 



Under these conditions, assuming Gaussian noise in each detector and optimal filtering, 
a filter function chosen to maximize the signal-to-noise ratio, SNR for two such detectors is 



7 We note that to allow for a small number of events that may be resolved at D C ~ 1, some authors 
conside r a more conserva tive threshold of DC ~ 10 to indicate the continuous regime (jRegimbau fc Hughes 



2009; 



Howell et al.ll201lh . 



8 We note however. IDrasco fc Flanagan! (|2003r ) also proposed a new statistic to search for SGWB with low 
duty cycles, which is currently being investigated in the LIGO/Virgo collaboration. 
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given by ( lAllen &: Romano! Il999l . equation 3.75) 



SNR 2 



9/7 4 

J -"0 rp 

50vr 4 



df 



/ 6 Pi(/)p 2 (/) 



(13) 



Here 7(/) is the 'overlap redu ction function', which accounts for the separation and relative 
orientation of the detectors (jFlanaganl 119931 ). and P%(f) and P%{f) are the noise power 
spectral densities of the detectors, and T is the integration time. As the optimal filter 
depends on Ogw(/)i a range of filter functions based on theoretical expectations of this 
function will need to be used. 

In this study we use data of relative positions and orientations for 10 i ndependent pairs 
of the five advanced detectors given in Table 3 of iNishizawa et al.l (120091 ) and employ the 
tensor-mode functions described in their equations (33-35). For ET we assume two detectors 
of triangular shape (60° between the two arms) and separate d by an angle of 120 °, for which 
the 7(/) has a constant value of —3/8 from 1 Hz to 1000 Hz (jHowell et al.ll201ll ). We adopt 
a value of SNR = 3 to indicate detection, corr esponding with false alarm rate of 10% and 
detection rate of 90% (lAllen &: Romano! Il999l ). We also assume an integration time of 3 
years for advanced detectors and 1 year for ET. 



5.3. Assessing the detectability using the worldwide network 

Calculating the SNRs for SGWB model (e) shown in Figure 3, we find a value of 
0.14 through cross-correlation by two the Advanced LIGO detectors H-L. For the other four 
models shown in Figure 3, we find variation in SNR of within 20%. For ET we find SNRs of 59 
and 112 assuming ET-B and ET-D sensitivities respectively for model (e), indicating that this 
signal will be easily-detected by third generation detectors. These results, based on average 
quantities, suggest that to detect the BBH background with two Advanced LIGO detectors 
will require a rate greater than even the higher rate estimate, r 2 ~ 0.43 Mpc~ 3 Myr _1 . As 
there will exist variation in the sensitivities, locations and orientations of detectors within 
a worldwide detector network, it is useful to compare the performances of different detector 
pairs and investigate how combining the network could improve the detection prospects. 

Two approaches of combining 2N dete ctors to increase t he sen sitivity of a stochastic 
background search have been proposed by lAllen fc Romano Jl999h . We apply these two 



methods to a network of 4 second generation detectors. In each case, the optimal SNR can 
be expressed as follows, with individual detectors (1-4) indicated in parenthesis: 

(i) Four-detector correlation (FC) - can be performed by directly correlating the outputs 
of 4 detectors 
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SNR 2 ptI « ( 12 )SNR 2 ( 34 )SNR 2 + ^SNR 2 ( 24 >SNR 2 + ( 14 )SNR 2 ^SNR 2 . (14) 

(ii) Combining multiple pairs (CP) - is performed by correlating the outputs of a pair 
of detectors, and then combining measurements from multiple detector pairs 

SNR 2 ptII = ( 12 >SNR 2 + ( 13 >SNR 2 + • • • + ^SNR 2 . (15) 

We now investigate the detectability of the BBH background by considering both these two 
approaches, as well as the cross-correlation method between two detectors. To do this, we 
substitute equation (TTU|) for f2 G w into equation (JT31 . Although there are uncertainties in 
ro and the true ranges of M c , this approximation allows us to quantify the advantages of 
the different approaches outlined above. The goal here is to provide insight into both the 
requirements for a detection and the constraints that can be supported by a null-detection. 

Firstly, we set M c = 8.7 M Q and compare standard cross-correlation measurements 
between different pairs of the five advanced detectors (H,L,V,C,A). We also determine any 
improvements that can be gained from the two approaches described above - (i) FC and (ii) 
CP. We assess performance through the values of r required to produce a SNR=3. 

Table [2] outlines the constraints on r$ we obtain through these different approaches for 
second generation detectors. We see that variation in the function 7(/) between different 
detector pairs together with different sensitivity levels, produces a variation in the values of 
r required for detection. Among the 10 independent pairs of the five advanced detectors, 
H-L performs best as indicated by the lowest required value of r (0.55 Mpc _3 Myr _1 ). We 
note that these values are above the higher rate r%. For cross-correlation measurements 
between two ETs, minimum detectable values of tq ~ 1.34 x 10~ 3 Mpc~ 3 Myr _1 and r$ ~ 
7 x 10~ 4 Mpc~ 3 Myr -1 are obtained using ET-B and ET-D sensitivities respectively. 

The values presented by Four-detector correlation, FC, suggest that no considerable 
advantage can be obtained through this approach. However, we see that by combining 
multiple pairs, CP, produces the most consistent improvement, the best of which comes 
from the AHLV combination (around a 40% improvement on H-L). Equation (15) shows that 
applying the CP method to the 3 identical arms of ET will reduce the minimum detectable 
values of r by a factor of a/3. 

5.4. The r —M c parameter space of a detectable BBH background 

We now adopt the same methods as employed in the last section to explore the parameter 
space of ro-M c required for a detectable BBH background. To allow for uncertainties in M c , 
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we consider a large range of (4-20) This range inc l udes m ost values of M c within the 



low metallicity (O.1Z ) distribution of iBelczynski et al.l (l2010bl ). We note however, that for 



the more realistic scenario that accounts for early common envelope mergers as the stars 
pass through the Hertzsprung gap, the range is more constrained, with M c around (4-9) 
M . 

Figure 5 shows the r -M c space for three different detection scenarios: 1) cross correla- 
tion of H-L - the best performing pair of advanced detectors; 2) AHLV (through CP) - the 
optimal combination of 4 advanced detectors; 3) cross correlation of two ET type detectors 
using two possible sensitivities (ET-B and ET-D). Zones encompassed by the two rates, r\ 
and r2, are shown by the shaded areas. 

The figure further confirms that advanced detectors are not likely to detect the BBH 
background at the rate r\. Detection would require th e high rate e stima te, r 2 and M c > 



1OM . Such a range for r -M c has been predicted by iBulik et al.l ( 120111 ) through studies 
of the two BH-WR systems, NGC300 X-l and IC10 X-l. A null detection would therefore 
confirm a lack of understanding in the binary evolution of such systems based on the single 
source models employed in this study. 

In regards to detection strategies for SGWB signals, Figure 5 shows some improvement 
in detectability through the use of CP. As already noted in Table [21 AHLV can improve the 
detectability by 40% compared with only two Advanced LIGO detectors (H-L). 

For ET, for the entire range of M c the signal would be comfortably detected at the lower 
rate r\. A detectable continuous background requires a rate of order 10~ 3 Mpc~ 3 Myr _1 or 
above. In comparison with ET-B, we find that a larger parameter space can be probed by 
ET-D due to an improved sensitivity at lower frequencies (< 20 Hz). 



6. Conclusions 



In this paper we have estimated the potential contribution of a population of coalescing 
BBHs to a SGWB sig nal. We are motivated by recent observations of BH-WR star systems 
(jCrowther et al.ll2010l) and by new estimates in the metallicity abundances of star forming 
galaxies ( jPanter et al.ll2008l ) that have suggested the rate of BBHs in field populations may 
be greater than previously expected. We base the single source emissions on energy spectra 
calculated from recent parameterized waveforms of lAjith et al.l (120081 . 120091 ) . Then, assuming 
that the BBH rate traces the SFR with some delay time, we derive cosmic source rate 
evolution models and extrapolate our single-source model out to high redshifts. Rather 
than a population synthesis approach, we assume average properties (e.g., masses and spins) 



-14- 



for the BBH population to determine the characteristics of the SGWB signal and principal 
parameters to which the BBH background is sensitive. 

Our results show that for M c < 1OM , the background is not likely to be detected 
through cross-correlation by two advanced detectors even at r 2 ~ 0.43 Mpc _3 Myr _1 . Only 
for greater values of M c , as have been predi cted to result from BH-WR systems such as 
NGC300 X-l and IC10 X-l ( TBulik et alJboilh . there is scope for detection. 



To further assess the detection prospects for second generation detectors, we have con- 
sidered the possibility of combining a worldwide network of advanced detectors to improve 
the cross-correlation statistic, namely the methods FC and CP described in section 5.3. 
We find that of these two approaches, CP can produce an improvement of up to 40% 
against a standard cross-correlation between two detectors. For the third generation de- 
tector, ET, the signal is accessible with a SNR of 59 and 112 at the lower rate estimate 
n ~ 3.1 x lO^Mpc^Myr" 1 using ET-B and ET-D sensitivities respectively. This signal 
could mask the primordial background signal at below around flaw ~ 4 x 10~ 10 at ~ 100 Hz. 

We note that the rates used in this study are computed assuming Milky Way type galax- 
ies and the standard formation channel - isolated binary evolution. Massive binary formation 
in early elliptical galaxies is expec ted to improve the coalescence rate (Ide Freitas Pacheco et al 
20061 ; lO'Shaughnessy et al.ll2010l ). This is particularly important for BBHs due to the longer 
delay time. In addition, dynamical formation scenarios in dense stellar environments can 
make a significant contribution to BBH rates. These other formation channels will not only 
increase the event rate of coalescing BBHs, but also add additional uncertai nty to the aver 



age co mponent masses of the BBH population. For example, simulations by lSadowski et al. 
( 120081 ) suggest that the average chirp mass in clusters is (M c ) ~ 2OM , much larger than the 
same found in the field (M c ) ~ 7M . This might indicate two similar backgrounds peaking 
at quite different frequencies. We have shown a minimum detectable r for ET at around 
10~ 3 Mpc ^Myr" 1 , which may be in the range of rate predictions from dynami cal formation 
scenarios ( jO'Leary et all 120071 ; ISadowski et al.l l2008t iMiller &: Lauburgi 120091 ). Therefore, 



detection at the sensitivity of ET could enable these two potential background signals to be 
untangled, thus allowing the average properties of the different populations to be probed. 

Additionally, clues to how these two populations contribute to a confusion background 
may be provided by Advanced LIGO/ Virgo through single detections in the shot noise 
(DC <C 0.1) regime. New data analysis techniques, such as the probability event horizon 
method ( jCoward fc Burmanll2005l ; iHowell et al.ll2007l ) which extracts the temporal signature 
from a population of transient sources or the maximum likelihood statistic (IDrasco fc Flanagan 
20031 ). could prove valuable in interrogating this regime. 
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In view of future searches for a primordial background signal, particularly for ET, sim- 
ulations incorporating both dynamical and isolated binary formation channels could prove 
useful. These will be considered in a future study. 
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Fig. 1. — Cosmic rate evolutio n models of BBH coa l escen ces. We show the effects of d i fferent 
SFRs by considering HB06 in Iflopkins fc Beacoml hood ) and W08 in Iwilkins et aD J2008h 
for a minimal delay time to = 100 Myr. To illustrate the effect of to we a l so show the case 
to = 500 Myr for the former SFR. 
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Table 2: Minimum values of the coalescence rate r (in Mpc _3 Myr _1 ) to detect the SGWB 
from coalescing BBHs assuming an average chirp mass M c = 8.7M . We explore the result 
from cross correlation measurements between different pairs of the worldwide network of 
second-generation ground-based detectors (H,L,V,C,A) and for two approaches of combining 
4 detectors (FC and CP), taking a SNR of 3 indicate detection. 
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Fig. 2. — GW energy spectra for a 10M Q -10M Q coalescing BBH in the non-spinning case 
and cases for non-precessing spins with three values of the single spin parameter x (see text). 
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Fig. 3. — The energy density parameter Oqw as a function of observed frequency for the 
BBH background corresponding to the five models given in Table 1. 
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Fig. 4. — The design sensitivity curves for future ground-based detectors: Advanced Virgo 
(Ad Virgo), LCGT, Advanced LIGO (Ad LIGO) and two possible configurations of the 
Einstein Telescope, ET-B and ET-D. 
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Fig. 5. — The detectable space of parameters for a stochastic background formed by coalesc- 
ing BBHs: the average chirp mass M c (in M & ) and the local rate density r (in Mpc~ 3 Myr -1 ). 
The r -M c curves correspond with those required to produce a SNR=3, through cross cor- 
relation measurements by two Advanced LIGO detectors (H-L) and two third generation 
detectors adopting ET-B /ET-D sensitivities, and the optimal combination of 4 advanced de- 
tectors (AHLV). The region above these curves can be considered the detectable parameter 
space. The two shaded re gions correspond t o zones encompassed by the two rates r 2 and ri, 
based on the estimates of iBelczynski et al.l ( l2010bl ). 



